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Abstract 

We present a quantitative description of the thermodynamics in a supercooled 
binary Lennard Jones hquid via the evaluation of the degeneracy of the inher- 
ent structures, i.e. of the number of potential energy basins in configuration 
space. We find that for supercooled states, the contribution of the inher- 
ent structures to the free energy of the liquid almost completely decouples 
from the vibrational contribution. An important byproduct of the presented 
analysis is the determination of the Kauzmann temperature for the studied 
system. The resulting quantitative picture of the thermodynamics of the in- 
herent structures offers new suggestions for the description of equilibrium and 
out-of-equilibrium slow-dynamics in liquids below the Mode-Coupling temper- 
ature. 
In recent years, a significant effort has been devoted to understanding the fundamen- 
tal nature of glass-forming materials, a long-standing open problem in condensed matter 
physics |I[]. Theoretical [§-0], experimental ^J^ and numerical efforts have broadened 
our knowledge of the physical mechanisms responsible for the dramatic slowing down of the 
dynamics in supercooled liquids (more than 15 order of magnitudes) as the temperature T 
changes over a small range, as the glass transition temperature is approached. 

Some recent theoretical approaches build upon ideas which were presented several 
decades ago [p|-[T^. In these works, the slowing down of the dynamics was connected to 
the presence of basins in configuration space. The short time dynamics (on a ps time scale) 
was related to the process of exploring a finite region of phase space around a local poten- 
tial energy minimum, while the long time dynamics was connected to the transition among 
different local potential energy minima [1^,1^. In this picture, upon cooling the intra-basin 



motion becomes more and more separated in time from the slow (and strongly T-dependent) 
inter-basins motion. The decrease of the entropy of supercooled liquids on cooling was 
associated with the progressive ordering of the system in configuration space, i.e. in the 
progressive population of basins with deeper energy but of lower degeneracy M. 



Following these ideas, Stillinger and Weber |12] introduced the concept of inherent struc- 
ture (IS), defined as local minimum configuration of the 3A^-dimensional potential energy 
surface. A basin in configuration space was defined as the set of points that — via a steepest 
descent path along the potential energy hypersurface — map to the same IS. This precise 
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operational definition of a basin allows the configuration space to be partitioned into an 
ensemble of distinct basins. Thus, the canonical partition function Zj^ for a system of A^ 
atoms at inverse temperature /? = l/ksT can be written as 

Zn = A-3^^exp(-/?<|.«) / exp[-/?A„(r^)]dr^ (1) 

where Ra is the set of points composing the basin a, $0 is the potential energy of minimum a 
and the non-negative quantity Aa(r^) measures the potential energy at a point r^ belonging 
to the basin a relative to the minimum. The integration over the momenta introduces the 



thermal wavelength A = Jph'^/27im, where m is the mass. Eq. (|l]) shows that both the IS 
energy and the thermal excitation within the basin region R^ contribute to Z^f. Stillinger and 
Weber also noted that, if the value of the potential energy minimum uniquely characterizes 
the properties of the basin, then a very powerful simplification of Eq. ([^) can be performed. 
By introducing a density of states ^{ejs) with IS-energy ejs, Zn can be written as 



^N 



deisn{ejs) exp[-(3eis - l3fi(3, ejs)] (2) 



where 



r dr^ 

PfiP, ers) = ln{ exp[-/3Ae,,$(r^)] — ^) (3) 

JRiers) A 



can be interpreted as free energy of the system when confined in one of the characteristic 
basins with IS energy ejs- Then, the probability that a configuration of the liquid extracted 
from an equilibrium ensemble of configurations at temperature T is associated to an IS with 
energy between ejs and ejs + Sejs is 

p, ^^ ^ ^{ejs)Sejsexp[-l3eis - PfjP, ejs)] ^ exp[-f3{eis - TSconfjeis) + f{(3, ejg))] 
^''''''^" ZNi(3) ~ ZNi(3) 

(4) 

where we have defined Sconfi^is) = kBln[Q{ejs)SGjs], since Q{ejs)S€.is is the number of 
states between ejs and ejs + Seis- 

The formalism proposed by Stillinger and Weber, although often used in the past to 
clarify structural issues in liquids [|T^,Q, for long time was not quantitatively applied to 
computer studies of the glass-transition problem, due to the significant computational effort 
required in equilibrating atomic configurations at low T. Only recently, Sastry et al fl^ 
addressed the problem of evaluating the T dependence of the average IS energy (e/g) in 
supercooled states in a binary mixture of L J particles, observing a significant decrease of eis 
on supercooling. This result, also observed for other models for liquids, e.g. in models for 
water [|l^ and orthoterphenyl |T^ , furnishes strong evidence of the relevant role played by 



the low-energy basins on cooling. In a recent work [^, we proposed to invert the relation 
between the e/5 energy and T to define an effective temperature at which the configurational 
part of the system is in equilibrium. This hypothesis, which has been proven useful in 
interpreting the aging process in a model liquid in terms of progressive thermalization of 
the IS [0], support the validity of Eq. (^ and together with the work of Sastry et al., calls 



for an effort in the direction of checking the formal expression for the supercooled liquid 
free energy, i.e. the T-range of validity, and an effort in the direction of evaluating the ejs 
dependence of the configurational entropy. This article is a first effort in this direction. 

The model system we study is the well-known 80-20 Lennard Jones A — B binary mixture 
ITSl , composed by a 1000 atoms in a volume Vo = (9.4)^, corresponding to a reduced density 
of 1.2039. Units of length and energy are defined by the a and e parameters of the A — A 
Lennard Jones interaction potential. The mass of atom A is chosen to be 1. In these 
units, ks = 1. Simulations, covering the range 0.446 < T < 5, have been performed in the 
canonical ensemble by coupling the system to a Nose'-Hoover thermostat ||T^. This system 
is well characterized and its slow dynamics has been studied extensively [|T^. The critical 
temperature of Mode Coupling Theory Tmct for this system is 0.435 [l^ . 



Between 500 and 1000 equilibrium configurations for each T (covering more than 80 
millions integration time steps for each T ) have been quenched to their local minima by using 
a standard conjugate-gradient minimization algorithm. By performing this large number of 
quenches we are able to determine not only ejs and its T dependence but also the probability 
distribution P{eis,T), shown in Fig.|l]-(A). 

In the T-region where Eq. (|) is supposed to hold, curves of ln[P{eis, T)] +I3eis are equal 
to Sconf{^is)lkB — Pf{P,eis), except for the T-dependent constant ln[Z{P)]. If f{(3,ejs) 
has only a weak dependence on e/5, then it is possible to superimpose P{ejs,T) curves at 
different temperatures which overlap in e/5. The resulting e/5-dependent curve is, except 
for an unknown constant, Sconf{eis)/kB in the e/5 range sampled within the studied T 
interval. This procedure is displayed in Fig.|l|-(B). We note that while below T = 0.8 curves 
for different T lie on the same master curve, above T = 0.8, curves for different T have 
different ejs dependence, thus showing the progressive e/5-dependence of f{P,eis)- The 
overlap between different P{ejs, T) curves below T = 0.8 indicates that the obtained master 
curve is indeed, except for an unknown constant, the e/5 configurational entropy, i.e. the 
logarithm of the number of basins with the same e/5 value. 

It is particularly relevant that, for T < 0.8, /(/?, e/5) ~ /(/3) and thus Z^ (Eq. |^) is well 
approximated by the product of a vibrational contribution [e~^'^] and of a configurational 
contribution depending only on the IS-energies and their degeneracy [J dejs^{,Gjs)G~^^'^]. 
Thus the liquid can be considered as composed by two independent subsystems, respectively 
described by the IS and by the vibrational part. The IS subsystem can be considered as 
a continuum of levels characterized by an energy value e/5 and an associated degeneracy 
Q{ejs)- When the IS-subsystem is in thermal equilibrium (with the vibrational subsystem), 
then the T dependence of the average configurational entropy Sconf can be evaluated using 
the standard thermodynamic relation 

dSconfiT) _ 1 

deis ~ T ^ ^ 

i.e. by integrating the T dependence of dejs/T. 

While the evaluation of the T-dependence of Sconf already furnishes relevant information, 
the evaluation of the unknown integration constant would allow for a determination of the 
number of IS with the same e/5 and, via a suitable low T extrapolation, to the determination 
of the so-called Kauzmann temperature T^, i.e. the T at which the configurational entropy 
appears to approach zero. To do so we exploit the fact that Eq. |] predicts that the liquid 
free energy FuquidiT) can be written as ||20[| 



Fu,u^d{T) = -kBTln[ZM{T)] = eis{T) - TS,onf{eis{T)) + /(/3, eis{T)) (6) 

where ejs{T) is the 67-5 value which maximises the integrand in Eq. (H). We assume that 
at the lowest studied T, the unknown /(/5, 67-5) can be approximated by the harmonic free 
energy of a disordered system characterized by the eigenfrequencies spectrum calculated 
from the distribution of IS at the corresponding T [^|. In this approximation, the en- 
tropy of the liquid, which can be calculate via thermodynamic integration, once the entropy 
of the corresponding harmonic disordered solid is subtracted, provides an estimate of the 
configurational entropy in absolute units. 

To calculate the liquid entropy we perform a thermodynamic integration first along 
the T = 5.0 isotherm, from infinite volume down to Vo, followed by a T integration of 
the specific heat at fixed volume, down to the lowest studied temperature. Since the T 
dependence of the potential energy E along the studied isochore is extremely well described 
by the law E{T) ~ T^/^, in agreement with recent theoretical predictions for dense fluids 
| |22| | , the potential energy contribution to the liquid entropy follows the law T~^/^. By 
summing the kinetic energy contribution we obtain an analytic expression which can reliably 
be extrapolated to temperatures lower than the studied ones. 

The evaluated T dependence of liquid and disordered-sohd entropies is reported in Fig.|^. 
We also show in Fig. ^ the contribution to the disordered-solid entropy arising from the 
T-dependence of the eigenfrequencies spectrum, to confirm that the T-dependence of the 
harmonic solid frequencies contribute only weakly to the entropy. The T dependence of 
AS = Siiquid-S disordered-solid is showu in Fig.^(A). We note that this difference vanishes at 
T = 0.297 ± 0.02, which defines Tk for the studied binary mixture. An independent recent 
estimate of T^ for this system, based on an integral equation approach and on a similar 
analysis of simulation data, is Tk = 0.29 p3|. The resulting ratio between Tk and Tmct 



support the view that the studied system has intermediate fragility character, as recently 
predicted by Angell and coworkers on the basis of a comparison between experimental results 



and numerical data for the same system [^ 



At the lowest studied temperatures, where the harmonic approximation is valid [^ 
AS" coincides with SconjiT). We use this identity to calculate the unknown integration 
constant for the inherent structure entropy. Moreover, by integrating Td{AS{T)), both the 
configurational energy dependence of the configurational entropy (Fig. |l|-(c)) and the T 
dependence of the configurational energy (0-(b)) can be evaluated, allowing to bridge the 
gap between Tk and the lowest T at which we were able to equilibrate the system. The 
present analysis predicts ejsiTK) = —7.82 ± 0.01. Thus, both the configuration entropy 
and energy around Tmct is halfway between Tk and the high T value, suggesting that the 
ordering process in configuration space at the lowest temperature which we have been able 
to equilibrate is far from being complete. 

The data reported in this article offer a quantitative thermodynamic analysis of the su- 
percooling state. This picture confirms the fruitful ideas put forward long time ago [^-p!^ 



and shows that a thermodynamic approach for the inherent structures subsystem becomes 
possible in supercooled states, since the thermodynamics of the inherent structures almost 
completely decouples from the "vibrational" thermodynamics. In particular, the quantita- 
tive evaluation of the degeneracy of the inherent structures for a well characterized system 
constitutes a basis for a comprehensive description of the slow dynamics below Tmct, the 
T-range in which an accurate theoretical prediction is still missing. Previous estimates of 



the degeneracy of e/5, the key ingredient for a descriptions of the dynamics in terms of IS, 
were available only for systems composed by less than 50 atoms [^,^. The proposed view 
of a supercooled liquids as composed by two weakly coupled subsystems — the IS subsystem 
and the "vibrational" subsystem — will also offer stimulating ideas for a microscopic un- 
derstanding both of the out-of-equilibrium thermodynamics theories recently proposed ||2^ 



and of the aging process [|l3)ill • 
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FIGURES 

FIG. 1. (A) Distributions P{eis,T) of the IS energy (per atom) for different equilib- 
rium temperatures T. From left to right: T = 0.446,0.466,0.5,0.55,0.6,0.8,1.0,2.0,4.0. (B) 
ln[P{eis,T)] + Pejs + C{T), for different equilibrium temperatures T. Symbols refer to different 
T-values. The curves which do not lay on the continuous line correspond to T = 4.0, 2.0, 1.5, 1.0, 0.8, 
from bottom to top. The constant C{T) has been chosen to maximize the overlap between curves 
with different T and the overlap with Sconfi^is) (in absolute units) calculated from AS{T), as 
discussed in the text. Sconfi^is) is reported as continuous line both in (B) and in the enlarged ejs 
range in (C). Note that when ejs = —7.82, Sconf = 0. 

FIG. 2. Disordered-solid (filled squares) and liquid (filled circles) entropies as a function of 
T. Filled rombs are S disordered- solid — 3A^/n(T/To), with Tq = 1, to demonstrate the weakness 
of the T contribution arising from the T-dependence of frequency distribution. Below T = 1, 
this weak T-dependent contribution is well fitted (long-dashed line) by the quadratic polynomial 
3201.62309 - 402.760885r + 199.228407T2, providing a rehable extrapolation of S disordered- solid to 
T lower than the studied ones (dotted line). The liquid entropy is extrapolated according to the 
theoretical predictions discussed in Ref. [^, i.e. Suquid + 1.5A;B/n(T/To) ~ T~^/^ 

FIG. 3. Configurational entropy (A) and energy (B) as a function of T calculated as difference 
of the liquid and disordered-solid entropies and energies, respectively. The arrow indicates Tmct 
for this system. The full circles in (B) are ejsiT). Note that in the region where the harmonic 



approximation for the disordered solid is expected to be valid |21|, eis{T) coincides with the 
configurational energy. 
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